The Poisson Rain Tessellation Under Spatial 
Expansion and Temporal Transformation 
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Abstract 



o 

For cell-division processes in a window, Cowan introduced four selection 
rules and two division rules each of which stands for one cell-division model. 
One of these is the area-weighted in-cell model. In this model, each cell is 
selected for division with a probability that corresponds to the ratio between 
its area and the area of the whole window. This selected cell is then divided 
by throwing a (uniformly distributed) point into the cell and drawing a line 
segment through the point under a random angle with the segment ending at 
the cell's boundary. 

For the STIT model which uses both a different selection and a different 
division rule, Martinez and Nagel showed that for a STIT process {Y(t, W) : 
t > 0} the process {e t Y(e t ,W) : t 6 M} is not only spatially stationary but 
also temporally. 

For a continuous-time area-weighted in-cell model it is shown by using the 
different and generalizing approach of a Poisson Rain that in order to get 



0^ ' temporal stationarity it is necessary to have an exponential spatial expansion. 

■ For the temporal transformation it is found that there is a strong relation 

between that transformation and the intensity of the Poisson Rain in time. 



MSC (2000): 60D05 

1 Introduction 

In PQ, Cowan examined cell-division processes in a window that are subdivided 
according to the ways a cell is chosen for division (by a selection rule) and then 
actually divided (obeying a division rule). Cowan considered area-, perimeter-, 
corner-number-weighted and equally-likely selection rules, where e.g. in the area- 
weighted case the probability for a cell to be chosen for division corresponds to the 
ratio between its area and the area of the whole window. The two division rules 
correspond to throwing a point into the cell or onto the cell's boundary (with the 
point uniformly distributed on the set where it is supposed to fall). It each case, a 
line is then drawn through that point under a random angle ending at the boundary 
of the cell. 

While Cowan examined his processes in discrete time it is not difficult to transform 
the discrete-time process into a continuous-time process by assigning a lifetime to 
each cell which is exponentially-distributed according to the weighting parameter of 
the cell (e.g. the lifetime of a cell C in the area- weighted case is £(A 2 (C))-distributed 
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with A2 describing the 

The STIT model which was introduced by Mecke, Nagel and Weifi and whose main 
characteristics are described in [5] fits into this categorization as an on-boundary 
Cowan model weighted according to the (generalized) perimeter. For a STIT process 
{Y(t, W) : t > 0}, Martinez and Nagel showed in [2] that the process {e t Y{e t , W) : 
t G M} is not only spatially stationary (as any STIT process) but also temporally. 
It is the goal of this paper to deduce necessary conditions for such a stability property 
to hold for the continuous-time area-weighted in-cell Cowan model. 
In order to do this, we first introduce a Poisson Rain with an intensity that may (or 
may not) vary over the time. It is quite straightforward to see that this Poisson Rain 
induces such a continuous-time area-weighted in-cell Cowan model as it does not 
matter whether we first choose a cell area-weightedly and then throw a uniformly- 
distributed point into the cell or whether we throw a uniformly-distributed point 
into the whole window and then consider the cell into which the point falls as the 
chosen cell. Of course, the lifetime distribution of a cell depends on the intensity of 
the Poisson Rain. 

To this end, we first examine the properties of the Poisson Rain and what happens 
when we expand the window in which it comes down and what happens when we 
transform the time. Then, conditions are deduced for the transformed Poisson 
Rain to have certain stability properties. Finally, it is explained how to come from 
the Poisson Rain to a tessellation; this tessellation in a sense inherits the stability 
conditions of the Poisson Rain and is called SAWSER if the stability conditions are 
fulfilled. 

It turns out that a SAWSER tessellation can only exist if the spatial expansion is 
exponential in time and if the Poisson Rain intensity and the time transformation 
fulfil a common condition. 

The idea for this paper came from a discussion with Richard Cowan, Werner Nagel 
and Viola Weifi on July 1, 2011. 



2 The Poisson Rain 

Let us first, independently of any tessellation, consider a Poisson point process. 
Let x : M 2 x [0, 00) — > [0, 00) be a differentiable function. Let further $ x be a 
marked Poisson point process on R 2 x [0, 00) x [0, tt) with an intensity measure 
x(x, u) ■ d(X 2 <E> Ai <8> V)(x, u, a) such that for any Borel subsets S C R 2 , T C [0, 00) 
and A C [0, tt) the number <& X (S x T x A) of points falling into S xT with an angle 
mark from A is determined by 

P($ X (S x T x A) = k) = 

F (/sxTxA^( z ' u )^ 2 ® Al ®'P){x,u,a)^ exp (- f SxTxA x( x , u)d(X 2 ® \\ ®V){x,u,a) 

(1) 

with Ai on the interval [0, 00) and A2 on M? being the one- and two-dimensional 
Lebesgue measure and V being a probability measure, respectively. 
Each point in $ x can be identified ) with Ti called the time of the point 

falling onto the plane and being the angle orthogonal to the line to be drawn 
through Xi. 

We will not need this angle in the next few sections but we will use it when we 
describe how we come to a tessellation from this point process. Therefore, we will 
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make use of the (unmarked) Poisson point process $ x which we can deduce from <3> x 
via the projection 

$ x (SxT) = $ x (5xTx[0,4 (2) 
The mean number of points of $ x falling into S during T is 

E$ x (SxT)=/ x(x,u)d(\ 2 ®\i)(x,u). 

JSxT 

This can also be written as 

E$ x (SxT)= f l SxT (x,u)x(x,u)d(\ 2 ®\ 1 )(x,u) (3) 

</R 2 x[0,oo) 

with Id being the indicator function for a set D. 

With % = 1, we get the usual homogeneous Poisson point process in space and time 
with intensity 1. 

We can interpret this Poisson point process as a Poisson Rain on the plane with 
the rain drops (the points Xj) touching the plane at the time Tj. The intensity of 
this rain in space and time is governed by the function %. In the following, we will 
consider the intensity of the rain as independent of the location, therefore we will 
write only x{ u )- Thus, we will now use the intensity measure x( u ) ' d{\2 <S> Ai)(x, it). 

3 Two approaches for expansion into the plane 

We now consider a convex and compact set 5cK 2 . We want this set to grow with 
the time and want to consider the mean number of points falling into the growing 
subset. Here, we have two possibilities to describe this growth until an arbitrary 
time t: 

• The first approach, to be called Growth During Rain (GDR), will see the 
set grow at the same time as the points are falling, so the growth of the set 
and the falling of points in the Poisson Rain are considered to take place 
simultaneously 

• The second approach, to be called Growth After Rain (GAR), will see first 
the points falling into a non-growing set until t and then grow the set to its 
size at time t. 

In both approaches, a point that has fallen into the set moves away from the origin 
as the time (and the set) grows. 

It should be clear that the GDR approach yields a stronger connection between the 
growth of the window and number of points falling into it than the GAR approach. 
Let us use monotonically increasing, differentiable functions £d : [0, oo) — > [0, oo) 
and £4 : [0, oo) — > [0, oo) for the GDR and GAR approach, respectively, to describe 
the growth of the set. Let us have a reference time Sq and Sq for which ^d(sq) = 1 
and ^a(sq) = 1 holds. Of course, this reference time may depend on the function 
£d or £a, respectively. Accordingly, by xd and xa we describe the intensities of the 
rain in the corresponding approach. 

Usually the time set T will now be an interval [s, s + 1] with s > and t > 0. 
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3.1 The GDR approach 

For the function £ D , let us define S = S s n = £ d (sq)S which is supposed to be a 
compact and convex set in the plane (a window). Thus we consider S to be our 
reference set with its size at the time Sq. Then we can define S s = £d(s)5 and 
S s+t = £d( s + t)S. For any time u G [s, s + 1], we can define S u = £d(u)S. Such a 
set has the area 



Also 



c c ( ,\ c £d(u) 



Thus, we can express the set S u (whose area grows with growing u) as a function 
of the set S s at the time s. A point (xj,Tj) of the Poisson Rain $ XD that hits the 
plane at the time hits the set S Ti if and only if X{ G S Ti . At the time s + t, this 
point has moved from xi to ^~z~y Thus, a point that fell at the time 73 is in S s+t 
if and only if it fell into S Ti . With the point process $ XD as defined for x through 
equations d2J) and ([T]) in hand, we can define for any time s with a set S^, a Borel 
set icl 2 and an interval [s, s + t) the point process 

which depends on the functions and £d, on the time s and the examined set S s 
as well as on the time s + t we consider. 

Thus, the mean number of points falling into that growing window S u = ^,d{u)S for 
u G [s, s + 1] can be calculated in analogy to (E} and considering the point process 
as a measure now as 

^ Xd ,Zn,s,S a ;s+t(S s+t X[s,3 + t]) 

= Jr2x[o,oo) 1 s u x[s, s +t](x,u)xD(u)d(X 2 <8) Ai)(x,m) 

= J[o,oo) (/ R a ls u (s)dA 2 (a;)) Xr»(w)dAi(u) 

= j S s +t ^2{S u )XD{u)du 

= f: +t ^X 2 (S s )x D (u)du 

= ^ i: +t e D (u)x D (u)du. 

All points of $x£>,?£>,s»&;s+t are P er definition contained within S s+t . 
This process is a Poisson point process: The points have moved away from the origin 
at the same speed as the window grows. This conserves the property that the mean 
number of points within a set is proportional to the area of the set. Additionally, 
all points have fallen independently of each other. 

3.2 The GAR approach 

Let us define again that S = S s a = £a(sq)S as a compact and convex set in the 
plane. Then we can say again that S s = Ca(s)S and S s+t = £a(s + t)S. 
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From the point process that is now denoted $ XA we can define a point process 

®XA,U,s,S;s+t = + t)Xi,Ti) : (Xi,Ti) G $ XA , Xi G S, 7} G [s,s + *]}. 

Note that the reference here is made to the set S which exists at the reference time 
s£ with U(so ) = 1- 
We then get 

E<f> XAtUjS: s; S+ t(S s+t x [s,s + t]) 
= E^ XA (Sx[s,s + t]) 

= X2 {iaJs) Ss ) r s +t XA{u)du 

= ir xA{u)du. 

The only difference to the GDR method is that the function i\(u) does not appear 
in the integral in contrast to i 2 D (u) in the GDR case. 

This process is a Poisson point process as well: The point process $ XA is a Poisson 
process which is only stretched to get the point process in question, thus keeping 
the Poisson properties of the unstretched process. 

4 Time transformation 

4.1 The function ip for both approaches 

We can now also transform the time. Let us have a function ip : (— oo, oo) — > [0, oo) 
which is differentiable on the support {set: ip{s) > 0} of ip. 
For the GDR method, by the £f, in the integral we have an intricate link between 
the number of points falling into the growing area and the area of the set ex- 
amined. Therefore, the only process that is a time-transformed GDR process is 
®XD£D,1>;8,s.,s+t which is defined by 

®Xn,Zn,M a , S +t(S s +t X [s,S + t]) = ^ XD ^ D ^(s),S^ s yMs+t)(^(s+t) X [lp(s), lp(s + t)]) 

for s G M and t > 0. We then get 

E ®XD£D,M„3+t(S a +t x [s,s + *]) 

For the GAR method, we have a fixed window S into which during the time interval 
[ip(s), ip(s + 1)] a certain number of points fall. Afterwards, we let the window grow; 
in this growth, the number of points in the growing set remains constant but the 
points themselves move away from the origin. Thus, we can set £ independent of ip 
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and define ® X A,U,il>;*,s,i>+t as 
^XA,^A,^;s,s,s+t 



{{£a(s + t)xi,Ti) : (xi,Ti) e ® XA ,Xi e S,n e [ifi(s),ifi(s + 1)]} 



®XA,UMs),S;f(s+t) 



(5) 



We then get 

^XA,U^\s,S,s+t(S s+t x [s,s + t)) 

= ^^XA 

,t A ,il>{s),S;iP{s+t)(S s+ t X [lfi(s),lfi(s + t)]) 

= HS) lZ +t) XA(u)du 

~ 15m his) X A \ u ) au - 
4.2 Equivalence 

If we start our consideration at an arbitrary time s (which is then transformed by 
ifi) and make the condition that the sets we are considering at that time (<SW S ) in 
the GDR, £a{s)S in the GAR approach as ifi does only have an influence on the 
number of points within the time interval but not the size of the set) have the same 
area, we can see from equations (0J and (JSJ) that with 

(&°0to = £to (6) 

for any s G R and 

(& ■ Xn)(u) = X a(u) (7) 

for any u G [0, oo) we have the same mean number of points falling into that set 
S$(s+t) for the GDR approach and S s+t for the GAR approach, respectively, until 
the time s + t. 

Note that from equation §6§ the requirement ^(sq ) = s o follows immediately. 



4.3 Conditions for the functions x-> £ an d ip 

Due to the equivalence of the approaches as per equations <$3§ and ([7]), in the fol- 
lowing, we will predominantly consider the GAR approach. In doing so, we will 
not always use the indices a and d when it is not necessary to make a distinction 
between the two approaches. 

In the first sections, we have stated that the functions x, £ and ifi shall be differen- 
tiable and non-negative. The requirement of differentiability will later make some 
calculations easier. Additionally, we want ifi to be monotonically increasing as it 
describes the time span we consider, thus ifi(s + 1) > ifi(s) should be true for t > 0. 
It should furthermore not be constant within the support {s G K : ifi(s) > 0} as 
then, every interval lying within such an area of constancy becomes a point yielding 
any integral zero. Thus, ifi is supposed to be strictly monotonically increasing on its 
support. Additionally, £ 2 should be monotonically increasing as otherwise speak- 
ing of a growing window would not make any sense; there also should be a time 
So G [0, oo) such that £(so) = 1- When transforming the time when we want also to 
have an s G (—00,00) such that ifi(s) = s and thus £(ifi(s)) = 1. 
This leads to the following conditions: 
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Condition 1 We only admit non-negative functions x> an d £ that are differen- 
tiate and for which 

^-i>{t) > for any t G R with i/j(t) > 

and 

holds. There also must be a time So G [0, oo) such that £(so) = 1 and a time 
s G (— oo, oo) such that i/;(s) = sq. 



5 Stability considerations 

For STIT tessellation processes, Martinez and Nagel showed in |2j that for a STIT 
tessellation process (Y t ,t > 0) the transformed process (e*y e *,t G R) is stationary. 
While we have not yet shown how a tessellation can arise from any of our two 
approaches, we can motivate two conditions for stability of our point process. 

• The first condition for stability is that the mean number of points per unit 
area becomes constant (and greater than zero) at least asymptotically. Other- 
wise, the process may degenerate by either having ever less or ever more points 
contained within a set or it may oscillate around a certain value without ac- 
tually approaching a limit. We will call that the MEPA condition below. 

• The second condition is that, given a set at an arbitrary time s, the waiting 
time for this (growing and moving) set to be hit does not depend on the time 
s. We will call this the INOT condition below. 

5.1 The MEPA condition 

The mean number of points per unit area (MEPA) is the ratio of the mean number 
of points in a set S s+t at a time s + 1 and the area of that set S s+t . Thus, we can 
write for some s G R 

- lim *^°o x^to 



llm ^°° &(8+t)\o(S„) ITUT U(s) XA[U)dU 



lim 



if that limit exists. This leads to 

i r^( s +*) 

I1L11L1UI1 £i iJtl lilt: blllblL lllil^oo TT 

ity 



Condition 2 Let the limit lim^oo ^ ^ +t ^ XA(u)du exist. Then, the inequal- 



t^oo £\(s + t) 

must hold for every s G R. 



< lim -p£- ( — J XA(u)du = M < oo 
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5.2 The INOT condition 



The independence of the observation time (INOT) condition is to ensure that when- 
ever we observe a (growing and moving) set the waiting time for a point of the point 
process to fall into it has the same distribution. 

Let us consider a set S s C M 2 at the time s. Then, the distribution of the waiting 
time Xs s for S s to be hit after the commencement (at the time s) of the observation 
can be described as 

(A (S ) 
-^TTT / XA ^ du 

(8) 

respectively. It is this waiting time's distribution which must be independent of the 
time s. Let us have a time Si and a set on one hand and another time s 2 and 
another set S^ on the other hand. Let the two sets have, however, equal area, i.e. 
A 2 (S«) = A 2 (S< 2 )). As A 2 (S«) = X 2 (U(si)S 1 ) = e A (si)HSi) for some set Si at 
the reference time s$ and, analogously, A 2 (S( 2 )) = A 2 (£a(s2)S 2 ) = £i(s 2 )A 2 (S 2 ) for 
some set S 2 at the reference time Sq, we can state for the sets Si and S 2 at the 
reference time 

A 2 (Si)^A 2 (S 2 )^ £4(^1) ^£4(^2). 

Because the sets and do have equal area however, we can state that the 
term 

W —- I XA(u)du 
must be independent of s. Therefrom, we can deduce 



Condition 3 The equation 



Ms+t) \ 
/ Xa{u)oIu = 



e A {s) 

must hold for every t > and every s6l. 

Note that as all functions are assumed differentiable the whole term is differentiable 
as well. 



6 The SAWSER tessellation 

We will first explain how we get a tessellation Y from the Poisson point process $ x . 
Then, the terms of a tessellation with an area-weighted selection rule (AWSER) and 
of a stable AWSER tessellation (SAWSER) will be defined. Lastly, we will prove 
that the tessellation Z arising from spatial expansion and temporal transformation 
of the original tessellation Y is indeed a SAWSER tessellation. 

6.1 The unexpanded and untransformed tessellation in W 

As the number of points falling into a bounded spatial set W during a bounded 
time interval [0,t] is almost surely finite, one can order these points according to 
the time they fall. Let k = & X (W x [0,t] x [0, 7r)) be the (random) number of 
points falling into W during [0,t] with any angle. Then we can with probability 1 
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re-arrange the indices in & X (W x [0,t] x [0,7r)) = {(xi, Ti, chi), (x K , r K , a K )} such 
that < ri < ... < t k < t. 

At the time Ti, the point xi falls into and a line segment 7i(xi, cci) D W is drawn 
through that point with the angle a\. Let us denote 

Y(r 1 ,W)= ll (x 1 ,a 1 )nW 

as the tessellation in W. Let us now define 

C(x 2 , y(n, w)) = cZ({y e w : xiy n y(r 1} = 0}) 

as the cell of Y(t 1: W) in which x 2 falls. Here, cl(S) means the closure of a set S 
and X2V is the connecting linear segment between x 2 and y. 
Then, we define 

y(t 2 , w) = F(n, w) u ( 7l (x 2 , a 2 ) n C(x 2 , y(n, w))) 

as the tessellation after the second point has fallen into the window. 
Thus, we iteratively define 

c(x n , y(r n _i, wo = cZ({3/ g iy : x^y n y(r n _i, = 0}) 

and 

y(r n , iy) = y(r n _i, u ( 7l (x„, «„) n c(x n , y(r n _i, iy))). 

Finally, under the given situation with k points in W during [0, t] we have 

Y(t,W) = Y(r K ,W). 



6.2 AWSER and SAWSER tessellations 

Let us now consider functions £ and tp as defined above. Then we define 

Z(t,W t )={(t)Y(m,W) 

as the tessellation arising from Y(ip(t), W) in a window W t = £(t)W and W t ) : 

t > 0) as the tessellation process. 

Let us now define in general the following. 

Definition 1 A tessellation Z = Z(t, W t ) in a compact and convex window W t = 
iit)W C M. 2 with non-empty interior is called AWSER if it employs an area-weighted 
rule for the selection of cells to be divided at each division step, i.e. if there is a 
cell division at a certain time t in the window W t , the probability for a certain cell 
C G Cells(Z(t,W t )) to be divided by the dividing segment S t is 

F(cnSt^a\w t ns t ^Hi) x - [l ' 



A 2 (W t ) 



Definition 2 An AWSER tessellation is called stable (and then a SAWSER tessel- 
lation) if 

• the ratio of the mean number of cells E\Cells(Z(t,W))\ and the area A 2 (Vy) 
of the window converges towards a positive constant for t — > oo: 

, E\Cells(Z(t,W t ))\ 
< lim -! ) )' t,n = M s < oo; 



9 



• for each cell C, its lifetime Xc is independent of the time of observation, i.e. 
if there are two cells C\ and C 2 of the same area (Xi{C-i) = X 2 (C 2 ) = L) at 
different observation times s% and S2, then 

F(X Cl > t\d G Cells(Z( Sl ,W Sl ), A 2 (d) = L) 

= F(X C2 > t\C 2 E Cells(Z(s 2 , W S2 ), \ 2 {C 2 ) = L) 
for any t > and any L > 0. 
Because of the Conditions [2] and El we can state 

Lemma 1 If the functions x> £ an d "0 fulfil the Conditions U\ on the regularity, [H 
(MEPA) and\M (INOT), then the tessellation Z{t,W t ) is SAWSER. 

Proof 

We have to prove that Z(t, W t ) is AWSER and that it is stable. Let us turn first to 
the AWSER property. 

For each point s in time, a cell C s has a lifetime Xq s for which the probability to 
exceed a time t can be described per equation (JSJ) by an exponential function with 
the cell's area at the observation time in the argument as 



-^7T / XA{u)du 

^AV S > Jib(a) 



Accordingly, the probability that the waiting time 

X Wa = min{X Cs : C s G Cells(Z(s,W s ))} 
of the state of the whole tessellation in W s to change exceeds t is 

P(X^ >t)= exp (-^T- ^(u)du 

\ t,A\ S ) Ji>{s) 

Due to well-known properties of the exponential distributions, the probability that 
this cell C s is the next to be divided is the ratio of its area and the area of the whole 
window W s . As the area of the cell and of the whole window grow proportionally, 
this ratio remains the same even within the growing window. This, however, means 
that the selection rule is area-weighted. 

The stability follows from the fact that for the number of cells \Cells(Z(t, W t )) \ the 
equation 

\Cells{Z(t,W t ))\ = <$> X {W x [0,ifi{t)) + 1 
holds. From Condition [2] we get 

,. nCeUs(Z(t,Wt))\ 

\i m E$ x (iyx[o#))+i 

- llUl^oo M{Wt) 

_ E<S> x (Wx[0,m) I U m 1 

— umt-j.oo \„(wa r nm^oo 



Aa(Wt) 1 """t^ 00 X 2 (W t [ 



02 



M + lim *->-^k) 



< 00 
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due to |£ 2 (s) > for any s > 0. 

Due to M > 0, of course also lim^oo E ^ Cell ^}^r t ^ Vt ^ > 0. 

A cell is a set in the window, so Condition [3] makes sure that the waiting time for a 
point to fall into that cell (i.e. the cell's lifetime) has a distribution which is inde- 
pendent of the time of observation. 

Thus, Z(t, W t ) is stable and has an area-weighted selection rule; it is a SAWSER 
tessellation. □ 



7 Dependencies between £ and ip 
7.1 A condition for £ 2 

We have the MEPA and the INOT condition in which we can re-write the integral 

as 

il>(s+t) 

XA (u)du = X(tP(s + t)) - X(^(s)) = (X o V0(s + 1) - (X o VOO) 

Ms) 

with X being the ant i- derivative of xa- 
Due to the INOT condition, we can say 

(Io^)( S + Mlo#)=/(^( S ) 

for some function /. As we must separate s and t, we need 

(Io^ + ^(Io^)( S )./» 

or Si0l = ]{t) for g = X o if) to hold. 
Differentiation to s leads to 

= — f(t) = d 9 ^ s + t "> = 9'{s + t)g{s)-g{s + t)g'(s) 
ds ds g(s) 9 2 ( s ) 

and thus 

g'(s + t)g(s)-g(s + t)g'(s)=0 
must hold for any t > 0. Then, however, 

g'(s) = g'(s + t) 
g(s) g(s + t) 

must hold for any t > 0. This means that 

— — = const. 

must hold. From the differential equation 

g'{s) = bg(s) 

we get 

g(s) = ae bs = (Xo^)(s) 
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for a, b G R. 

So, we get the equation 

ae Ks+t) _ Qe 6, = /(f)£( s ) 

and thus because of the condition £a(so) = 1 

f{t) = ae bso (e bt - l) and $(s) = e - bs °e bs . 
Putting this into the MEPA condition, we get that 

< lim e - b{s+t) a (e b{s+t) - e bs ) = lim a (l - e~ bt ) = a = M < oo. 

Thus, the MEPA condition is fulfilled as well for b > 0. 

It is noteworthy that b = 0, i.e. = 1 for all t > 0, is not admissible as the MEPA 

condition would not hold then (the limit would be zero as (Xoif))(s+t) — (Xo^))(s) = 
0). 



7.2 Conditions for \ an d ip 

Per the condition 

(Xoif])(s) = ae bs 

we get 

X (^))^W = a6e ta - 



As ^V ; (' s ) > 0) we can define ^ 1 as the inverse function of if> such that for any s 

ds ' 



with ^(s) > we have ip 1 (ip(s)) = s. We will write -^ip(s) = if>'{s). Then we have 



lf]'( S ) 

and then, inserting s = 

X{t)=x{ ^ m= ^- Wy (9) 

This is e.g. fulfilled for (Case A) 

^( s ) = ae bs and x(s) = 1 (10) 

(due to ip" 1 ^) = \ In ^) or for (Case B) 

V>(s) = sl[ ,oo)(s) and x(s) = afre fe 

due to , 0~ 1 (s) = s = if>(s) on the support of if>. Here, 1 is the indicator function 
which is necessary to have if) defined on whole K. but map into [0, oo) only; such a 
if] is differentiable on its support as required. 

For the classification Cowan introduced in [TJ, the special case xi s ) = 1 corresponds 
to the model Cowan-2a, namely the area-weighted model with in-cell division rule. 
As we can deduce from equation ffTUj) . we must transform the time by if>(s) = ae bs 
in order to get a stable tessellation. 
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7.3 Theorem 

Theorem 1 A tessellation 

z(t,w) = at)Y(m,w) 

with Y(t, W) arising from a Poisson point process $ x is SAWSER if and only if the 
conditions 

£ 2 (t) = e - bso e bt 

for b > and £(sq) = 1 and 

x ^(s))^(s) = abe bs 
as 

are fulfilled. 

For any given ip fulfilling the Conditions, 

abe^' 1 ^ 

x(t) 



is the fitting rain intensity. 

7.4 The re-translation to the GDR process 

For the re-translation to the GDR process, we can use equations flO]) and (J7J). We 
have a ip which only must be strictly monotonically increasing, £ A (t) = e - bso e bt and 
X as described per equation (j^J). 
We can then quite straightforwardly see that 

and 

, . ab 

XD(t) 



Let us consider two combinations of functions (x,£ 2 sVO; namely (l,e , e ) and 
(abe bt , e bt , tl[o t00 )(t)) , i.e. functions with so = 0. We then get 

(Case A) ( Xa , ^, = (MV 6t ) ^ ( %D ,^, = (ar 1 ,*, e") 

and 

(CaseB) = (abe bt , e bt , tl [0>oa) (t)) <^ ( X d, 4>){t) = {ab, e bt , tl [Q>c 

It should be noted that in Case A in the GDR approach x is n °t defined for t — 0. 
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